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Abstract 

This paper proves that in Size Theory the comparison of multidimensional size functions 
can be reduced to the 1-dimensional case by a suitable change of variables. Indeed, we show 
that a foliation in half-planes can be given, such that the restriction of a multidimensional 
size function to each of these half-planes turns out to be a classical size function in two 
scalar variables. This leads to the definition of a new distance between multidimensional 
size functions, and to the proof of their stability with respect to that distance. 



Keywords: Multidimensional Size Function, Multidimensional Measuring Function, Natural 
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Introduction 

Shape comparison is probably one of the most challenging issues in Computer Vision and Pat- 
tern Recognition. In recent years many papers have been devoted to this subject and new 
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mathematical techniques have been developed to deal with this problem. In the early 90's, Size 
Theory was proposed as a geometrical/topological approach to shape comparison. The main 
idea is to translate the comparison of two datasets (e.g. 3D-models, images or sounds) into the 
comparison of two suitable topological spaces M and N, endowed with two continuous functions 
ip : M — > R fc , ip : N ^ R fe _ These functions are called k- dimensional measuring functions 
and are chosen according to the application. In other words, they can be seen as descriptors 
of the properties considered relevant for the comparison. In |19| the definition of the natural 
pseudo-distance d between the pairs (M, (p), was introduced, setting d (jJA, ip), (3\f, ip)j 

equal to the infimum of the change of the measuring function, induced by composition with 
all the homeomorfisms from M to 3\f. Unfortunately, the study of d is quite difficult, even 
for k = 1, although strong properties can be proved in this case (cf. [§1 Hip. Size Theory 
tackles this problem by introducing some mathematical tools that allow us to easily obtain 
lower bounds for d, such as size homotopy groups and size functions (cf. |19| and |10j1- The 
idea is to study the pairs (M,(p ■< x),JA(<p ^ y)), where M,((p ■< t ) is defined by setting 
M(<p <t) = {P e M : (pi{P) < U,i = l,...,fc} for t = (ti,...,t k ) € R k . The fc-th size 
homotopy group ilk (x, y) describes the non-trivial equivalence classes of fc-dimensional loops in 
JA(tp ^ x) that remain homotopically non-trivial also in M.(ip ^ y). Size functions count the 
number of connected components in M.(tp ^ y) that meet M.(tp ^ x). It turns out that ttq(x, y) 
is a set whose cardinality is equal to the value taken by the size function at (x,y). From the 
homological point of view, an analogous approach, named Size Functor, has been developed in 
|2] for 1-dimensional measuring functions. 

More recently, similar ideas have independently led Edelsbrunner et al. to the definition of 
Persistent Homology (cf. |12l I13j). and Allili et al. to the definition of the Morse Homology 
Descriptor (cf. £Q). 

From the beginning of the 90's, size functions have been studied and applied in the case 
of 1-dimensional measuring functions (cf., e.g., [HI [SI [THI CHI [IHl I2DI EH 1221 El)- The 
multidimensional case presented more severe difficulties, since a concise, complete and stable 
description of multidimensional size functions was not available before this work. 

In [3], Carlsson and Zomorodian examine the completeness problem by studying Multidi- 
mensional Persistent Homology. In that paper, it is claimed that multidimensional persistence 
has an essentially different character from its 1-dimensional version. Indeed, their approach 
does not seem to lead to a concise, complete and stable descriptor in the multidimensional case, 
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whereas it does in classical Persistent Homology (see 0]). 

The first result of this paper is the proof that in Size Theory the comparison of multidimen- 
sional size functions can be reduced to the 1-dimensional case by a suitable change of variables 
(Theorem^. The key idea is to show that a foliation in half-planes can be given, such that 
the restriction of a multidimensional size function to these half-planes turns out to be a clas- 
sical size function in two scalar variables. Our approach implies that each size function, with 
respect to a fc-dimensional measuring function, can be completely and compactly described by 
a parameterized family of discrete descriptors (Remark |2J. This follows from the results proved 
in |17| about the representation of classical size functions by means of formal series of points 
and lines, applied to each plane in our foliation. An important consequence is that we can 
easily prove the stability of this new descriptor (and hence of the corresponding /c-dimensional 
size function) also in higher dimensions (Proposition 0), by using a recent result of stability 
proved for 1-dimensional size functions (cf. 013) and analogous to the result obtained in 0] 
for Persistence Homology. As a final contribution, we show that a matching distance between 
size functions, with respect to measuring functions taking values in K fe , can easily be introduced 
(Definitional). This matching distance provides a lower bound for the natural pseudo-distance, 
also in the multidimensional case (Theorem GJ . All these results open the way to the use of 
Multidimensional Size Theory in real applications. 

Outline. In Section ^ we give the definition of fc-dimensional size function. In Section [5] the 
foliation we use is presented, and the reduction to the 1-dimensional case is proved. Section [21 
shows the stability of our computational method, implying a lower bound for the natural pseudo- 
distance. Additionally, a new distance between multidimensional size functions is introduced. 
In Section^ the effectiveness of the multidimensional approach is tested on an example. Section 

examines some links between multidimensional size functions and the concept of vineyard, 
recently introduced in 0. Sectional concludes the paper, presenting some discussion and future 
work. 

1 Definition of k-dimensional size function 

For the present paper, M, N denote two non-empty compact and locally connected Hausdorff 
spaces. 

In Multidimensional Size Theory any pair (M, (p), where <p = ((pi, . . . , ipk) ■ M — > K fc is a 
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continuous function, is called a size pair. The function tp is called a k- dimensional measuring 
function. The following relations ;< and -< arc defined in R k : for x — {x\, . . . ,x k ) and y — 
(yi,...,y k ), we say x -< y (resp. x -< y) if and only if x t < y { (resp. x t < y { ) for every 
index i = 1, . . . , k. Moreover, R fc is endowed with the usual max-norm: \\(x\, X2, ■ ■ ■ , x k )\\oo = 
maxi<i<fc \xi\. In this framework, if M and N are homeomorphic, the size pairs (M.,0) and 
(N, tp) can be compared by means of the natural pseudo-distance d, defined as 

d((M,0,CN,$) - inf max||^(P) - $/(P))||oo , 

where / varies among all the homeomorphisms between M and N. The term pseudo-distance 
means that d can vanish even if the size pairs do not coincide. Here, and in what follows, R fe x R fe 
and R 2k are identified. 

Now we introduce the /c-dimensional analogue of size function for a size pair (M, 0). We 
shall use the following notations: A+ will be the open set {(x , y) € R k x R fe : x -< y}, while 
A = <9A+. For every fc-tuple x = (xi, ...,x k )e K fc , let M(<p d x) be the set {P e M : <fii(P) < 
Xi , i — 1 , . . . , ky . 

Definition 1. For every k-tuple y = (j/i, . . . , y k ) € R k , we shall say that two points P,QeM 
are (0 ^ y) -connected if and only if a connected subset ofM(0 ^ y) exists, containing P and 
Q. 

Definition 2. We shall call (fc-dimensional) size function associated with the size pair (M, 0) 
the function t(M,tp) '■ A + ~~ * N, defined by setting £(M,ip)(x, y) equal to the number of equivalence 
classes in which the set M(0 < x) is divided by the (0 -< y) -connectedness relation. 

Remark 1. In other words, £(M,<p)(x,y) counts the connected components in M(0 ^ y) con- 
taining at least one point ofM{0< x). 

2 Reduction to the 1-dimensional case 

In this section, we will show that there exists a parameterized family of half-planes in R fe x R fe 
such that the restriction of £{M,<p) to each of these planes can be seen as a particular I-dimensional 
size function. 

Definition 3. For every unit vector I = (l\, . . . ,l k ) ofR k such that h > for every i = 1, . . . , k, 
and for every vector b = (61, . . . , b k ) of R k such that Y^i=i h — 0, we shall say that the pair 
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(I, b) is admissible. We shall denote the set of all admissible pairs in R fc x R fc by Adrrik- Given 
an admissible pair (l,b), we define the half-plane TiVff) of R k x R k by the following parametric 
equations: 

x = si + b 
y = tl + b 
for s,t6R, with s < t. 

Proposition 1. For every (x,y) G A+ there exists one and only one admissible pair (l,b) such 
that (x, y) G 7T (; - g) . 

Proof. The claim immediately follows by taking, for i = 1, . . . , k, 

, yi-Xi _ xi J2j=i Vi - y^ Ej=i x o 



Then, x = si + b, y = tl + b, with 



£*=i*j Ej=i(%-^) 



Now we can prove the reduction to the 1-dimensional case. 



□ 



Theorem 1. Let (I, b) be an admissible pair, and Ff rj \ : M — ► R be defined by setting 

{1,0) 



FZ~(P)= max 



¥>i(P)-&i 



' i=l,..,fc [ 

T/ien, /or every (x, y) = (si + b,tl + b) G ""(f g) ^ e following equality holds: 

t(M,0){x,y) = ^(M,F^ S) )( S '*) • 

Proof. For every x = (ari, . . . , xt) G M fe , with = sZ, + b{, i = 1, . . . , k, it holds that M(<p ^ 
i*) = JVUi 7 '^- < s). This is true because 

(Lb) 

M(^f) = {PeM:<pi(P)<Xi,i = l,...,k} = 

= {P e M : tpi(P) < sh + b u i = l,...,k} = 

= {p 6M :^<M = l t}=M(^<.). 



Analogously, for every y — (t/i, . . . , y^) € R fc , with yi = tli + bi, i = 1, . . . , fc, it holds that 
3Vt(<p <y) = M(f€ ? < t). Therefore Remark [T] implies the claim. □ 

In the following, we shall use the symbol F%j\ m the sense of Theorem ^ 

Remark 2. Theorem^ allows us to represent each multidimensional size function as a param- 
eterized family of formal series of points and lines, on the basis of the description introduced in 
\17j for the 1-dimensional case. Indeed, we can associate a formal series with each ad- 

missible pair (/, 6), with C(fg) describing the 1-dimensional size function F <p y The family 

{°(fS) ' C'^) ^ ^4d?7ifc| is a new complete descriptor for l(M,tp)> i n the sense that two multi- 
dimensional size functions coincide if and only if the corresponding parameterized families of 
formal series coincide. 



3 Lower bounds for the k-dimensional natural pseudo- 
distance 

In d , it has been shown that 1-dimensional size functions can be compared by means of 
a distance, called matching distance. This distance is based on the observation that each 
1-dimensional size function is the sum of characteristic functions of triangles. The match- 
ing distance is computed by finding an optimal matching between the sets of triangles de- 
scribing two size functions. For a formal definition we refer to E] (see also 0] for the 
analogue of the matching distance in Persistent Homology). In the sequel, we shall denote 
by d ma tch{£cM p <? \A, 3 J the matching distance between the 1-dimensional size functions 

The following result is an immediate consequence of Theorem ^ an d Remark 

Corollary 1. Let us consider the size pairs (M,(p), (N,^). Then, the identity ^(m,^) = ^(h$) 
holds if and only if d ma tchUcM w „i J = 0, for every admissible pair (l,b). 

The next result proves that small enough changes in (p with respect to the max-norm induce 
small changes of I p ^ . with respect to the matching distance. 

Proposition 2. If (M, if), (M, \) are s ^ ze pairs and maxp e jvt ||<£>(P) — x(P)||oo < e? then for 
each admissible pair (l,b), it holds that 

e 
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Proof. From the Matching Stability Theorem 25 in [5] (see also 0) we obtain that 

Let us now fix P S JVC. Then, denoting by t the index for which max^ y '^^ b ' is attained, by 
the definition of Ff-~. and we have that 

(Lb) (l,b) 



P% (P) - F* (P) = max -^-^ - max 

^ t -(P) - Xi(P)-bi , n{P)-h xt(P) - h 
T t m f X k " T t T t 

iptjP) - xtjP) K MP) - x(P)l|co 



min i= i ... fc k 



In the same way, we obtain P*L (P) - F*L (P) < ilggh^-^Jlis . Therefore, if max PeM \\0{P) - 
x{P)\\oo<e, 

\0(P)-X(P)\\oo s 



max 

PeM 



F^fP)-F* (P) 

(l,b) y ' (l,b) y ' 



< max ■ 



psm min,=i,...,fe k mint=i,...,k h 

□ 



Remark 3. Analogously, it is easy to show that small enough changes in (l,b) with respect to 
the max-norm induce small changes of £ (M with respect to the matching distance. 

Proposition |3 and Remark [3] prove the stability of our computational approach. 
Now we are able to prove our next result, showing that a lower bound exists for the multi- 
dimensional natural pseudo-distance. 

Theorem 2. Let (JVC, <p) and (N, ip) be two size pairs, with M, Ji homeomorphic. Setting 
d((M, ip), (JSf, tj))) = inf/ maxpgjvt \W{P) ~ V'(/(-P))llooj where f varies among all the homeomor- 
phisms between M and N 7 it holds that 

sup min k-d match (£ (MF )J m ^ J < <p), (3\f, V))). 

^6)eAdm fc i - 1 '-- !fc ' < ! - b » k ' 

Proof. For any homeomorphism f between M and N, it holds that £, ,? 

Moreover, by applying Proposition[3with e = maxpgjvt ||<?(P) ~~ ^(/(-P))l|oo an d X = 4* fi an( i 
observing that o / = F%°J , we have 

(;,6) J (l,b) 



i=i 



min /i • d ma tch(e {MtF 0_ y £ mF t J < max ||£(P) - 



for every admissible (l,b). Since this is true for each homeomorphism / between M and J\f, the 
claim immediately follows. □ 
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Remark 4. We observe that the left side of the inequality in Theorem defines a distance 
between multidimensional size functions associated with homeomorphic spaces. When the spaces 
are not assumed to be homeomorphic, it still verifies all the properties of a distance, except for 
the fact that it may take the value +00. In other words, it defines an extended distance. 

Definition 4. Let (M, (p) and (N, tp) be two size pairs. We shall call multidimensional matching 
distance the extended distance defined by setting 

(r,S)eAdm k l - 1 '---' k ' ('• b ' 1 ' (W 

Remark 5. If we choose a non-empty subset A C Admk and we substitute sup^- jj- jGAdnik with 
sup,fj\£ A in Definition^ we obtain an (extended) pseudo-distance between multidimensional 
size functions. If A is finite, this pseudo-distance appears to be particularly suitable for applica- 
tions, from a computational point of view. 



4 An example 

In R 3 consider the set Q = [— 1, 1] x [— 1, 1] x [— 1, 1] and the sphere § of equation x 2 +y 2 + z 2 = 1. 
Let also $ = ($1, $2) : M 3 ^ R 2 be the continuous function, defined as $(x, y, z) — (|ar|, \z\). In 
this setting, consider the size pairs (M, (p) and (N, tp) where M = 9Q, N = §, and <p and ijj are 
respectively the restrictions of $ to M and 3\T. In order to compare the size functions £(m,0) an( i 
a, we are interested in studying the foliation in half-planes i>fgs> where I = (cos 9, sin 9) 
with 9 £ (0, §), and b — (a, —a) with set. Any such half- plane is represented by 



< 





= s COS t 


t + a 


%2 


= s sin 6 


- a 


V\ 


= t cos 9 


+ a 


!J2 


= t sin 9 


— a 



with s,t € R, s < t. Figure n shows the size functions £,~, „<? . and I x , for 9 = j 
and a = 0, i.e. / = and 6 = (0,0). With this choice, we obtain that i^jL = 

V / 2max{( y 5i, £^2} = \/2max{|a;|, \z\} and = v / 2max{V'i, ^2} — V2max{\%\, \z\}. Therefore, 
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Theorem n implies that, for every (xi, X2, J/i, 1/2) G ^ffg- 



(i,6) 

s s i i 
s s i t 



V (Lb)' 



In this case, by Theorem [3] and Remark[5] (applied for A containing just the admissible pair that 
we have chosen), a lower bound for the natural pseudo-distance d((M.,<p), (N, ijj)) is given by 



matchiii-M w <? J is equal to the cost of moving the 



Indeed, the matching distance d 
point of coordinates (0, y/2) onto the point of coordinates (0, 1), computed with respect to the 
max-norm. The points (0, V2) and (0, 1) are representative of the characteristic triangles of the 



In other words, the 



size functions £^ M F j and £^ ^ y respectively. 

We conclude by observing that £(M,<pi) = ^(H,i/>i) an d £(M,<p 2 ) = ^Wtfa) 
multidimensional size functions, with respect to <p, ip, are able to discriminate the cube and the 
sphere, while both the 1-dimensional size functions, with respect to tfi,(p2 and ip\,ip2, cannot 
do that. The higher discriminatory power of multidimensional size functions motivates their 
definition and use. 



z 




' y 


/ 


-1 




















I 


X 








/ 













Figure 1: The topological spaces M and N and the size functions I 



sin 



■ 



CH,F% _ ) 



associated 



with the half-plane 7t,,-a , for I = ( 



V2 V2 



) and b = (0,0). 
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5 Links between dimension vineyards and multidimen- 
sional size functions 



In a recent paper jS], Cohen-Steiner et al. have introduced the concept of vineyard, that is a 
1-parameter family of persistence diagrams associated with the homotopy ft, interpolating be- 
tween /o and /i . These authors assume that the topological space is homeomorphic to the body 
of a simplicial complex, and that the measuring functions are tame. We shall do the same in 
this section. We recall that dimension p persistence diagrams are a concise representation of the 
function rank H^ v , where H^ v denotes the dimension p persistent homology group computed 
at point (x, y) (cf. Therefore, the information described by vineyards is equivalent to the 
knowledge of the function rank Hp ,y , computed with respect to the function f t . We are inter- 
ested in the case p — 0. Since, by definition, for x < y, rank Hq' v coincides with the value taken 
by the size function in^j t ^{x,y), it follows that, for x < y, dimension vineyards contain the 
same information as the 1-parameter family of size functions {^(M,/ t )}t6[o,i]- Anyway, another 
interesting link exists between dimension vineyards and multidimensional size functions. This 
link is expressed by the following theorem. In order to prove it, we need the next two lemmas. 
The former states that the relation of (0 < y) -connectedness passes to the limit. 

Lemma 1. Assume that (M, 0) is a size pair and y = (yi, . . . , y^) € R fc . //, for every e > 0, P 
and Q are (0 ^ {y\ +£,..., yk + s) ) -connected in M, then they are also (0 ^ y) -connected. 

Proof. For every positive integer number n, let K n be the connected component of M(</? -< 
("Hi +)•••) Vk + )) containing P and Q. Since connected components are closed sets and M 
is compact, each K n is compact. The set |~) K n is the intersection of a family of connected 
compact Hausdorff subspaces with the property that K n+ \ C K n for every n, and hence it is 
connected (cf. Theorem 28.2 in P- 203). Moreover, f] n K n is a subset of M.(tp ^ y) and 
contains both P and Q. Therefore, P and Q arc (0 ^ y )-connected. □ 

The following lemma allows us to study the behavior of multidimensional size functions 
near A (where they have not been defined because of instability problems when the measuring 
functions are not assumed to be tame). 

Lemma 2. Let (M, 0) be a size pair. If ' x -< y then lim e ^ + ((m,^) i^, {y\ + e, . . . , yk + s)) is 
equal to the number L(x,y) of equivalence classes of M(p x) quotiented with respect to the 
(0 ^ y) -connectedness relation. 
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Note that, for x -< y, L(x,y) simply coincides with £(m,<p) (^iV)- 

Proof of Lemma\Q First of all we observe that the function i(M.,ip) (jJi + e i ■ ■ ■ > Vk + e )) IS 
nonincreasing in the variable e, and hence the value lim £ ^ + £(M,ip) (j/i + £ i ■ ■ ■ > Vk + e)) is 
defined. The statement of the theorem is trivial if lim e ^ + £(m,<p) (vi + e i ■ ■ ■ > + e )) = 
+00, since, for every s > 0, the inequality £(m,<p) (2/1 + e i ■ • ■ i Vk + e)) < £(^7 y) holds by 
definition, and hence the equality L(x,y) = +00 immediately follows. Let us now assume that 
lim e _ f0 + i(M,<p) (Vi + £ i ■ ■ ■ 1 + e )) = r < +oo- hi this case a finite set T = {Pi, . . . , P r } of 
points in M(<^ d a?) exists such that, for every small enough e > 0, every P € JA((p < x) is 
(¥> d (Vi +£,•■• , y/b + e))-connected to a point Pj G 3 in M. Furthermore, for i ^ j the points 
Pi, Pj G y are not ^ (j/i +£>•••> S/fc + e))-connected and hence not (tp d y )-connected either. 
From LemmaHit follows that every P € M(</? ^ af) is (0 < y)) -connected to a point Pj G 3 in 
M. Therefore y) = r = lim e ^ + ^(3YW) (yi + E > • • • > Vk + £))• □ 

Theorem 3. Por i G / = [0, 1], consider the family of size pairs (M, ft) where f t is a homotopy 
between f : M -> R and /1 : M -> R. De/me x : M x J -> R 3 6y x(P, i) = (ft(P), t, -t). Then, 
for every t € I and x,y el with x < y, it holds that 

rankH^ v (t) = lim i{Mxi,xj(x, i, — i, y + e, i + e, — i + e), 

where Hq' v (t) denotes the dimension persistent homology group computed at point (x,y) with 
respect to fj. 

Proof. We know that rank Hq' v (t) is equal to the number of equivalence classes of M(/( < x) 
quotiented with respect to the (ft < y) -connectedness relation. On the other hand, Lemma [3 
states that lim e ^ + £(Mxi,x) ^ — ^> 2/ + £ j ^ + e > + e ) is equal to the number of equivalence 
classes of M x I(x r< (8, f, — £)) quotiented, with respect to the (x di (y,i, —t)) -connectedness 
relation. By definition of x, this last number equals the number of equivalence classes of M(/f < 
x) quotiented, with respect to the (ft < y)-connectedness relation. This concludes our proof. □ 

However, although these two links exist, the concept of multidimensional size function has 
quite different purposes than that of vineyard. First of all, vineyards are based on a 1-parameter 
parallel foliation of R 3 , while the study of multidimensional size functions depends on a (2k — 2)- 
parameter non-parallel foliation of A + C R fe x R k . In fact, multidimensional size functions 
are associated with /c-dimensional measuring functions, instead of with a homotopy between 
1-dimensional measuring functions. Furthermore, |Hj does not aim to identify distances for 
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the comparison of vineyards, while we are interested in quantitative methods for comparing 
multidimensional size functions. 

6 Conclusions and future work 

In this paper we have proved that the theory of multidimensional size functions can be reduced 
to the 1-dimensional case by a suitable change of variables. This equivalence implies that multi- 
dimensional size functions are stable, with respect to the new distance D matc h, and this allows 
us to use them in concrete applications by exploiting the existing computational techniques. 
Many theoretical problems deserve further investigation, among them we list a few here. 

• Choice of the foliation. Other foliations, different from the one we propose are possible. 
In general, we can choose a family T of continuous curves 7*3 : M — > R k such that (i) for 
s < t, 7«j(s) ~< 7a(£), (ii) for every (x, y) G A+ there is one and only one 75 G T through 
x, y and (iii) the curve 75 depends continuously on the parameter a (this last hypothesis 
is important in computation for stability reasons). It would be interesting to study the 
dependence of our results on the choice of the foliation. 

• Extension to the algebraic context. We think that the main results obtained in this 
paper for multidimensional size functions can be straightforwardly extended to the ranks 
of size homotopy groups and persistent homology groups. 
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